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ABSTRACT 

Symbiotic X-ray binaries (SyXBs) comprise a rare class of low-mass X-ray binaries. 
We study the Galactic SyXBs, which we consider as detached binaries composed of 
low-mass giants and wind-fed neutron star companions, by simulation of the interac- 
tion of a magnetized neutron star (NS) with its environment and utilizing a popula- 
tion synthesis code. We focus mainly on the parameters that influence observational 
appearance of the SyXB: the donor wind velocity (vw) and the angular momentum 
distribution in the shell of matter settling onto NS. We estimate the birthrate of SyXB 
as ~ 4.1 x 10~ 5 yr _1 to ~ 6.6 x 10~ 6 yr" 1 and their number in the Galaxy as ^(100 
- 1000). Assumed stellar wind velocity from cool giants is the input parameter that 
influences the model SyXBs population most. 

Among known SyXBs or candidate systems, 4U 1954+31 and IGR J16358-4724 
in which NS have very long spin periods may host quasi-spherically accreting NSs. 
GX 1+4 has a peculiar long-term spin behaviour and it may also be a quasi-spherical 
wind- accreting source. We cannot identify whether there are wind-fed accretion disks 
in 4U 1700+24, Set X-l, IRXS J180431. 1-273932 and 2XMM J174016.0-290337. 
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. . ■ 1 INTRODUCTION 

> ' 

. ,—i ' Usually, symbiotic stars are binaries in which a hot white 
dwarf (WD) component burns hydrogen accreted from a 
5-h cool giant companion which loses mass via stella r wind 
.Cy.' or Roche lobe overflow (|Tutukov fc Yungelson] 1 197fJ ). The- 
oretical studies of the population of symbiotic stars have 
been published, e.g., by lHan et all (fl995l ); Yungelson et al.l 
( 19951 ): llben fc Tutukovl l|l99rJ ): lLii et al.l (|2006l . 120081 . 120091 



12011J). 

Recently, a small but rapidly growing subclass of symbi- 
otic stars gained attention, in which accreting component is 
a neu tron star (NS ) . All these systems are hard X-ray emit- 
ters (|Miirset et al.| [l997). They are lo w-mass X-ray bina - 
ries, with possible exception of Set X-l jKaplan et al ]|2007l) 
for wh ich a high-mass solution also is viable. iMasetti et al.l 
(2006) dubbed these systems symbiotic X-ray binaries 
(SyXBs). Currently, 10 binaries are classified as SyXBs or 
candidate systems (Table [TJ. 

Remarkable features of SyXBs are their orbital periods, 
among the long est determined for low-mass X-ray binaries 
|Liu et all [20071 ) and long NS spin periods. With a 18300 s 
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spin period, 4U 19 54+31 is the slowes t accretion-powered 
NS known to date (Mas etti et al.l l2006). Long spin periods 
of NS in SyXBs are not easy to understand assuming the 
standard disk accretion. In that case, equilibrium NS spin 
periods would require magnetar-like magnetic fields of NSs. 
However, long periods find natural explanations if accretion 
onto NS proceeds quasi-spherically. The quasi-spherical ac- 
cretion onto NS wa s reconsidered in the recent paper by 
IShakura et all (|2012l ). It was shown that at small X-ray lu- 
minosities, L x < 4 x 10 36 erg/s, the subsonic settling accre- 
tion regime from the stellar wind onto NS takes place, dur- 
ing which an effective removal of angular momentum from 
the NS magnetosphere can take place through hot convec- 
tive shell formed above the magnetosphere. In this accre- 
tion regime, equilibrium spin periods of NS can be very 
long even for the standard values of the NS magnetic fields 
~ 10 12 - 10 13 G. This theory was applied for explanation 
of the observed features of low-luminosity X-ray pulsars, in- 
cluding long-term spi n-down of NS and spin- luminosity cor- 
relations in GX 1+4 (Gonzalcz-Galan c talll2012l ). the 5-h 
spin period of 4U 1954+31 dMarcu et al.ll201ll). the NS spi n 
and luminosity behaviour in X Per ( Lutovinov et al.ll2012l ). 
etc. 

Clearly, the model of the Galactic population of SyXBs 
needs a detailed treatment of accretion process onto NSs, to 
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Table 1. Parameters of observed symbiotic X-ray binaries. Columns 1 to 7 list the name of the star, spin period P s and its deriva- 
tive P/P s , orbital per iod P or b, X-ray luminosity, the distance from the Sun, spectr al type of companion to NS. References: C97- 
IChakrabartv fc Roche! il997h; B9 7- lBildsten et al.Nl997l); H06 - lHinkle et al.l 12009) ; G12-fGonzalez-Galan et all |2012l '); C08- ICorbet et al. 
teOOST ); M06 -lMasetti et all d200^ ; M02-lMasetti et al.H2002l) ; ~K07jKaplan et al.1 1120071 ); M07-lMasetti et , al.1 hoOW ; N07jNucita et al. 
1 20071') ; P04 jPatel et alj i200^~P0 7- |patel et alj j20Q7h; B0 6- lBodaghee et alj fe006lh~ N10 jNespoli et alj j2010t l; L05- lLutovinov et al. 
1 20051 ); T06 jThompson et alj i2006t) : FlO- lFarrell et alj ll20ld) ; Mll jMasetti et alj <201ll) . 
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which we pay particular attention in the present paper. We 
focus on the population synthesis for SyXBs. We simulate 
zero-age population of NSs with companions and then follow 
evolution of the population with special accent on the spin 
and magnetic field of NSs interacting with the matter lost 
by their companions. In §2 we present our assumptions and 
describe some details of the modelling algorithm. In §3 we 
present a detailed example of the evolution of a NS inter- 
acting with its companion. In §4 the properties of the model 
population of SyXBs are presented and individual observed 
SyXBs are discussed. Conclusions follow in §5. 



2 THE MODELS 

For the simulation of bin ary evolution, we us e rapid binary 
st ar evolution code BSE (IHurlev et al.ll2002T ) with updates 
bv lKiel fc HurlevH 20061 ). If any input parameter is not spe- 
cially mentioned, its default value is taken from these papers. 
The code is extended with an algorithm for the treatment 
of spin evolution of a magnetized NS, as described in §J09 
|2~3~T1 [2321 

Like in the main case considered in our study of s ymbi- 
otic stars with white dwa rf components l)Lii et al.ll2006l ). we 
use iMiller fc Scalol l|l979t) initial mass-function for primary 
comp one nts, a flat distribution of initial mass ratios of com - 
ponents dKraicheva et all 1 19891 ; iGoldberg fc Mazehl [1993) . 
We assume that all binaries have initially circular orbits. 
After a supernova explosion, new parameters of the orbit 
are d erived using standard formulae, (e. g., IHurlev et al.l 
12002th 

The model is normaliz ed to formation of on e binary 
with Mi > O.8M per year ijYungelson et al.lll993h . We use 
10 7 binary systems in the Monte-Carlo simulations. This 
gives a statistical error ^5% for the number of Galactic 
SyXBs. 

Below, we only mention several specific assumptions 
used in the code. 



2.1 Kick Velocity 

Nascent NS receive an additional velocity ( "kick" ) due to a 
still enigmatic process that violates spherical symmetry dur- 
ing the collapse of a massive star or later. The kicks have 
dichot omous nature, as it was sugges ted quite early bvlKata 
ll 19751) and l ater c onfirmed by, e. g.. lHartman et all (| 1997t l : 
IPfahl et all lj2002t h Observationally, the kick is not well 
constrained due to numerous selection effects. Currently, 
high (~ 100 km s _1 ) kicks are associated with NS origi- 
nating from core-collapse supernovae, while low kicks (~ 
10km s _1 ) — with NS born in electron-capt ure supernovae 
and a ccretion-induced collapses (see, e. g., Ilvanova et al.l 
(2008) for refe rences and discussion ). 

Following Podsiadlo wski et al.1 (|2004T ) . we assume that 
core collapses are experienced by stars with ZAMS mass 
M/Mq ^ 11. For NSs born via core-collapse we apply 
Maxwellian distribution of kick velocity i>k 

= ,/^4e-" 2/2 < (1) 

V 7rcr£ 

We use velocity dispersion <7k = 190 km s~ , which 
is consistent with the data on pulsar proper motions 
( Han sen fc Phinnevl[T997l ). Variation of <7k between 50 and 
200 km s _1 , introduces an uncertainty < 3 in t he birthrate 
of low- and intermediate-mass X-ray binaries (|Pfahl et al.l 
2003). Since in this paper we focus on the physical param- 
eters of close binaries that mostly affect observational ap- 
pearance of NS — their equilibrium spin period and X-ray 
luminosity, we do not discuss the effects of <Tk on SyXB's 
population. 

Electron-capture supernovae (ECS) and accretion- 
induced collapses (AIC) are associated with formation of 
ONeMg cores of stars or white dwarfs. Coll apse is trig- 
gered by electron captures on 20 Ne and 24 Mg (|Mivaii et al.l 
1980). ECS occur in single stars, while AIC in white dwarfs 
in close binaries. The range of stars that form ONeMg 
cores or dwarfs depends on subtle details of the treat- 
ment of rotation, mass loss, mixing etc. and also depends 
on stellar evolution code applied (se e, e. g., di s cussio n in 
iPodsiadlowski et al.|[20o3 ). Following Kiel et al.1 (|2008l ). we 
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assume that ECS occurs if stellar helium core mass is 1.4 ^ 
Mhb/Mq ^2.5. Progenitors of these stars have ZAMS mass 
between ~ 8.0 and 11.0 M Q ijHurlev et alJl2000l ). Accretion- 
induced collapses happen when accreting ONeMg WD reach 
the Chandrasekhar limit. Progenitors o f ONeMg WDs hav e 
ZAMS mass between ~ 6.3 and 8.0 M e (|Hurlev et al.ll2000h . 

Con sidering the proble m of retention of NSs in globular 
clusters. IPfahl et alj 1120021 ) suggested that NSs born in ECS 
and AIC have velocity dispersion lower than 50 km s _1 . In 
this study we assume it equal to 20 km s . 

Response of ONeMg dwarfs to accretion is t reated in 
the s ame way as for CO WD (see for the details iLii et all 
12009ft . 



2.2 Mass Transfer and Angular Momentum Loss 

Mass transfer may occur by accretion of the matter lost via 
stellar wind or due to Roche-lobe overflo w. If the system is 
detached, we apply iBondi fc Hovkj (1 19441 ) accretion formula 
in which mass-accretion rate greatly depends on the wind 
velocity i> w . D etermination of u w is difficult. According to 
lHarperi ||l996ft . the main characteristic of the cool winds of 
evolved K and early M giants is that their terminal velocities 
are lower than the surface escape velocity, typically, « w ^ 
l/2v esc . As the wind velocity is very important parameter 
for accretion, we set w w = l/2v eS c and u w = 2i> csc in different 
simulations. 

For the case of RLOF, the critical issue is whether mass 
loss occurs on dynamical time scale and common envelope 
(CE) forms. This depe nds on the mass rat io of components 
q = M d /M a . Following fHurlev et all (|2002l ). we assume that 
condition for formation of CE (if the accretor is a main- 
sequence star) is q > 4 for the donors in the Hertzsprung 
gap. If the donor is a giant, critical condition is given by 
q > q C r, where 



<?cr = 0.362 + ■ 



1 



3 x (1 - M c /M) ' ^ 

wher e M c is the mass of the stellar degenerate core l|Webbinkl 
1 19881 ). 



For common envelopes, we apply IWebbinkl (f l984) equa- 
tion for energy balance, as modified by|deKoo| (jl990|) by 
adding a numerical "structure" factor A, meant to describe 
the dependence of binding energy of the donor on the den- 
sity distribution. Then, final separation of the components 
after CE-stage depends on the product of A and "common 
envelope efficiency" a cc (the fraction of binary binding en- 
ergy which is spent to expell com mon envelope). We assume 
Q CC x A=0.5. IPfahl et all (|2003h have shown that such high 
value of a ce x A is necessary for explanation of the birthrate 
of Galactic population of LMXB. Reduction of q cc x A may 
result in decrease of the birthrate of SyXB due to increase 
of the number of mergers. 

In the treatment of angular momentum loss we follow 
original BSE code wit h only except i on: fo r magnetic braking 
we replace Eq. (50) oflHurlev et al.l |2002l ) by a formula from 
IVerbunt fc ZwaarJ (fl98ll ) with A = 1 and k 2 = O.fl 



The reason for this change is the following. iHurlev et all i2002l ) 
use a mo dification of magnetic b raking law parametrization sug- 
gested bv lRappaport et al.l lll983l ). Then, BSE produces a long- 



2.3 Evolution of Neutron Stars 



The regime of interaction of a rotating magnetized NS 
(single or in a binary system) with surrounding mat- 
ter has been recognized to be the most i mportant astro- 
hysical aspect of its evolution, se e , e. g., [P ringlc fc Reea 
1972h ; llllarionov fc Sun^aevl (ll975l); lGhosh fc Lambl (|l978l . 



ll979j |bh~ lLovelace et al l j 19951 . 1 19991 )). From the point of 
view of population synthesis, a conv enient way of d escrib - 
ing NS evolution was elaborated by iLipunov et all fl992). 
Progress in the 3D MHD simulations of the interaction of 
a rotating magnetized N S with accreting plasma generally 
confi rmed the basic ideas (|Romanova et al.ll20o3 .[2003. 2004, 

lioosl) . 

The regime of interaction of a rotating magnetized NS 
with its environment can be determined by relations between 
four characteristic radii: 

• the radius of the light cylinder R\ = c/w where c is the 
speed of light and u) is the spin frequency of NS; 

• the radius of the gravitational capture (the Bondi ra- 
dius) 



Rg = 



2GM NS 



where Mns is the mass of NS and Woo is the velocity of 
surrounding matter; 

• the corotation radius R c = (GM/cj 2 ) 1/3 ; 

• the radius _R s t where the flow of accreted matter (free- 
fall or accretion disk) is stopped due to interaction with the 
NS magnetosphere. 

Depending on the relations between the above four 
radii, NS can be in several basic evolutionary states 
l|Lipunov et al.lll992l) : 

• Ejector state [Rat > max(i?G, Rl)]> in which the pres- 
sure of relativistic particle wind is sufficient to keep the 
stellar wind plasma of the companion away from the NS 
magnetosphere. For example, ordinary radio pulsars form a 
subclass of ejectors (additional physical conditions may be 
required for effective relativistic plasma creation); 

• Propeller state [J? c < R st < max(i?G, Ri)], in which the 
relativistic wind pressure cannot prevent the infalling mat- 
ter from interaction with the magnetosphere, while the fast 
rotating magnetic field prevents the matter from stationary 
accretion onto the NS surface; 

• Accretor state (i? a t < Rc), in which accretion onto mag- 
netized NS is centrifugally allowed, and the NS magnetic 
field channels the accreting material toward the magnetic 
poles of the NS. 

We shall assume that i? st is equal to the Alfven radius 
Ra- For disk and quasi-spherical supersonic accretion (see 
below) 



Ra = (m 2 /(^nsV2GMns)) 2/7 , 



(3) 



where Mns is the accretion rate onto NS, magnetic dipole 
moment fi = Bns-Rns/2, and Bns and _Rns are the magnetic 
field and the radius of NS, respectively. 

However, Eq. (|3]) is valid only if the accretion rate is 



living population of systems with giant donors with mass 
O.lMg, contradicting observations. 
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not too high and the accretion luminosity is below the Ed- 
dington limit. For t he su percritical accretion via a disk, 
IShakura fc Sunvaevl l|l973r i suggested that the accretion rate 
in the disk begins to fall monotonically from certain radius 
R — R s (the spherization radius). The super-Eddington 
accretion onto magnetized NSs was studied i n det ail by 
iLipunovl (|l982T l and summarized in ILipunovl (|l987t ) ■ The 
spherization radius _R S is approximately given by 



Rs 



4-7TC 



Mns, 



where n is the opacity of the accreting matter. When Mns 
is higher than the Eddington accretion rate, NSs can accrete 
via the disk at a rate 



M NS = 



Ra 
R* 



M NS . 



The excess of the matter is expelled in the form of a wind 
from a supercritical accretion disk and carries away the spe- 
cific orbital angular momentum of the NS. In the s upercrit- 
ical re gime the Alfven radius is determined as (see ILipunovl 



Rsa = 



V I~<a'2(Vm) 



fJ, 2 K 



2/9 



(4) 



Similarly, in the case of super-Eddington accretion, depend- 
ing on the relations between R\, Rg, Rc, and Ra, NS can 
be in super-ejector ( SE), super-pro peller (SP) and super- 
accretor (SA) states (jLipunovll 19871 ). 

For a NS, "evolution" means that its basic parame- 
ters which determine the interaction with the surrounding 
medium change. First of all, this concerns NS spin and mag- 
netic field. The effect of the NS mass increase in the accre- 
tion state is generally less important, unless a hypercritical 
accretion (especially inside common envelopes) is allowed. 
Then the accret ion-induced fo r matio n of a black hole may 
occur (see, e. g.. lBrown fc Led l|2004 ) for more details). 



2.3.1 Spin evolution 

Spin evolution (spin-up or spin-down) of a NS in a binary 
system can be conveniently described by an angular momen- 
tum conservation equation 



= K BU — K B d , 



(5) 



where I is NS momentum of inertia, K au and K s< ± are spin- 
up and spin-down torques, respectively. Their forms are dif- 
ferent in different evolutionary states and also depend on 
the mode of accretion, i.e., whether accretion disk is present 
or accretion proceeds quasi-spherically. At small accretion 
rates, spin evolution of a NS in a binary system proceeds 
almost identically to that of a single NS, so we neglect in- 
teraction with matter for M < 10~ 15 M o yr" 1 and treat NS 
spin evolution as in the ejector state. 

In the ejector state, K au = and NS spins down due 
to current braking iBeskin et al.l[l9 93). To within a numer- 
ical factor depending on the angle £ between the spin and 
magnetic dipole axes, the spin-down torque can be written, 
using the light-cylinder radius R\, as K s a = 2/i 2 /(3i? 3 ). The 
secular change of the angle £ is beyond the scope of our 
consideration. 



In the propeller state, accretion is centrifugally pro- 
hibited, K su = 0, and the spin-down torque can be writ- 
ten using the magnetospheric radius in the form K B a = 
fct/i 2 / i?A- The numerical f actor, k t is model-dependent but 
is ~ l jShakura et alj|2012h . 

In the accretor state, both K su and K s a are determined 
by the geometry of the accretion flow, as specified below. 

Roche lobe overflow. If accretion onto the compact star 
occurs via Roche lobe overflow through the vicinity of the 
inner Lagrangian point, an accretion disk is formed 0. The 
spin-up torque in both accretor and super-accretor states is 
determined by the specific angular momentum at the inner 
disk radius and is K BU — Mns V GMnsRa • The spin-down 
torque for disk accretion in the first approximation can be 
written in all cases using the corotation radius R c as K B a = 
(l/3)^ 2 /i? c 3 . 

The total torque exerted on a NS in the accretor or 
super-accretor state is then 



K = M ns VGM ns Ra - l/3/x 2 7? c - 



(6) 



The competitive action of spin-up and spin-down torques, 
on average, diminishes the total torque acting on NS, and 
an equilibrium stat e of NS is reached. The eq uilibrium NS 
spin period is then l|Lipunov fc Postnovll 19881 ') 



P s oq = 5.72M, 



-5/7 



M 



-3/7 6/7 



(7) 



'NS J "16 ^30 

where /i3o = /i/(10 30 G cm 3 ) is the NS dipole magnetic mo- 
mentum, and Mi6 = Mns/(10 g/s) is the accretion rate. 

For super-Eddington accretion, the NS equilibrium pe- 
riod is found to be 



1.76 x 10"V3o 3 ^n S 2/3 



(8) 



Wind-fed accretion. In the wind-fed X-ray binaries, the 
physical condition for formation of an accretion disk is j a > 
j'k(Pa), where j a = k v Q^,R G is specific angular momentum 
of the captured stellar wind matter, f2b = 27r/P or b is orbital 
frequ ency, fc w is a numerica l coeffi cient, which we set to 0.25 
after llllarionov fc Sunvaevl l|l975h . Jk(Ra) = VGMnsRa is 
the Keplerian angular momentum at the NS magnetosphere. 
If an accretion disk forms, the case of wind- fed disk accretion 
is realized, otherwise, a regime of quasi-spherical accretion 
is established. We shall assume that the wind-fed disk ac- 
cretion is similar to the disk accret i on via Roche lobe. 

As discussed bv lShakura et all |2012l ). there can be two 
different regimes of quasi-spherical accretion. First, if the 
gas heated up in the bow shock cools down rapidly, the mat- 
ter falls toward NS supersonically. This happens if accretion 
rate is fairly high (above ~4x 10 16 g/s). In this regime, all 
gravitationally captured matter eventually reaches the NS 
surface. Spin evolution of NS is determined by the sign of 
the specific angular momentum of captured matter, which 
can be positive (prograde) or negative (retrograde), and NS 
can sp in-up or spin-d own, respectively. Numerical simula- 
tions (iHo et al.l Il989h show that the sign of the specific 
angular momentum of the captured matter can alternate; 
here, however, we shall consider only the prograde case, 



2 Unless the mass-donating component turns out to be inside the 
magnetosphere of the compact object; this case can be relevant 
for the accretion onto a magnetized white dwarf, as in AM Her 
systems 
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i.e. only NS spin-up during the supersonic accretion with 
K su = MnsVGMRa- To within a factor of 2, the magneto- 
spheric radius Ra in this regime coincides with that for disk 
accretion [Eq.([3|]. and here we shall ignore the difference 
(but see below). 

Second, if accretion rate onto NS is smaller than 4 x 
10 16 g/s, shocked gas remains hot, a quasi-static extended 
shell forms around the NS magnetosphere, and accretion 
proceeds in the subsonic settling regime. In this regime, the 
actual accretion rate onto NS Mns can be lower than the 
capture rate by NS, and its value is essentially determined 
by the ability of plasma to enter the NS magnetosphere via 
instabilities; however, in this paper we shall assume that in 
this regime all matter captured by the NS from the stellar 
wind ultimately accretes onto the NS, as in the case of su- 
personic accretion. This assumption implies that the X-ray 
luminosity at the model SyXB stage can be higher than in 
reality, le ading to shorte r equi librium NS spin periods. As 
shown by Sha kura et al.l l|2012f ). large-scale convective mo- 
tions set in the shell, and the latter mediates angular mo- 
mentum transfer to or from the NS magnetosphere depend- 
ing on X-ray luminosity. The NS spin-up/spin-down equa- 
tion can be written in the form given by Eq.©, with the 
spin-down torque 



K ed = Z2Mnsw.Ra 
and the spin-up torque 



(9) 



(10) 



Here, the dimensionless coefficient Z\ determines the 
plasma-magnetosphere coupling and it is a function of the 
accretion rate Mns and NS dipole magnetic moment /x, 
Z2 — Zi — 2/3. The value of the index n reflects dif- 
ferent rotational distributions of matter inside the shell 
(oj{R) ~ R~ n )\ it depends on the treatment of (generally, 
anisotropic) turbulent viscosity responsible for angular mo- 
mentum trans port in the shell. Th e analysis of observed 
X-ray pulsars ( Sh akura et al.l I2012T ) suggests that an iso- 
angular momentum distribution with n = 2 is most likely, 
but quasi-Keplerian law with n — 3/2 cannot be excluded 
as well. Therefore, in our modeling we shall consider both 
quasi-Keplerian (n = 3/2) and iso-angular momentum dis- 
tributions (n = 2) as representative cases. 

It is important to note that the value of the Alfven 
radius itself in the regime of settling accretion is different 
from that in the case of disk accretion or supersonic quasi- 
spherical accretion (Eq. Q above): 



Ra w 10 h 



j4o\ 2/n ( M NS \ 
MieJ \1-5M@J 



-2/11 



(11) 



After inserting this expression into the formulae for spin- 
up /spin-down torqu es and using the value of Z\ from 
IShakura et all (|2012r ), the NS spin evolution equation reads: 



Iu 



■ 3+2n 

AM— 



BM 



3/11 



(12) 



The accretion-rate independent coefficients A and B are (in 
CGS units): 



.4 : 
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(13) 



B « 5.4 x 1(T 
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(14) 



Here Vs = Wu,/(1000 km/s) is the relative stellar wind ve- 
locity, K\ is a dimensionless numerical factor which differs 
for different systems; we shall assume K\ = 40, which is 
the typi cal value found from the analysis of observed X-ray 
pulsars (|Shakura et al.|[2012T ). 

In the equilibrium state (when the total torque exerted 
on the NS vanishes, K au + K B d — 0), the equilibrium spin 
period of a NS in the settling accretion state is independent 
of the source-specific coefficient K\ and it is given by 

P^s3 1000- (0.034r- 2 4^i n Mr 6 * 

/P \ ( M \ 13/H(2-n)-26/ll (15) 
x ( p n,h ] / Mns ] „ 

In our calculations, we take Mns = 1.5M© in Eqs. (11), 
(13)-(15). Among other factors, this one appears to be 
known quite well. In all formulae we omit some lengthy fac- 
tors describing anisotropic turbulence in the quasi-spherical 
envelope; taking them into account would change the result 
to within less than 10 per cent. 

We also stress that in the free-fall supersonic accretion 
case which is realized for high X-ray luminosities (L x > 
4 x 10 36 erg/s), spin evolution of a NS is entirely deter- 
mined by the sign of the captured specific angular momen- 
tum, and, strictly speaking, the equilibrium period cannot 
be determined from the mass accretion rate and NS mag- 
netic field only; in real systems the equilibrium state still 
could be established due to alternating spin-up/spin-down 
torques in the inhomogeneous stellar wind, but this cannot 
be taken into account in our simulations. For prograde angu- 
lar momentum accretion, the NS in this regime will spin-up 
until the corotation radius will become comparable with the 
magnetospheric one, Ra = Rc, most likely resulting in the 
transient accretion or propeller regime. Transient sources are 
beyond the scope of the present paper. 

The values of spin-up and spin-down torques exerted 
on a rotating magnetized NS are summarized in Table [3] 
The spin evolution of the NS with these torques at each 
state can be explicitly understood as a tendency to reach an 
equilibrium period, P eq , which is derived from the condition 
K — 0. Equilibrium periods are also listed in Table [2] 

As mentioned in £12.21 mass loss and and mass exchange 
can change the orbital period. In the ejector and propeller 
states, the matter transferred from the companion is ejected 
out of binaries. In this work, the ejection of matter is treated 
as occurring in Jeans mode, i.e. as spherical wind ejection 
with the specific angular momentum equal to the NS orbital 
one. 

Since overwhelming majority of NS in SyXB form via 
core-collapses, we set initial spin periods of neutron stars 
equal to 10 ms l|Cordes et al.ll2004 ). 



2.3.2 Magnetic field evolution 

For the decay of the magnetic field of an accreting 
NS, a commonly accepted idea d oes not exist as yet. 
iBisnovatvi-Kogan fc Komberd (|l974f ) suggested that the de- 
cay results from the scre ening of the origin al magnetic field 
by the accreted matter. iRudermar] (| 19911 ) suggested that 
the decay is due to the crustal motion on the surface of NS. 
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Table 2. Spin-up K BU and spin-down K sd torques acting on the rotating magnetized NS at basic evolutionary states, and the correspond- 
ing equilibrium periods derived from the equation K BU — K sd = 0. We use units fi^o = n1tn ^ :i , M\p = Mns 
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1000 ■ (0.034) r 



M N3 
1.5M Q 
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iGoldreich fc Reiseneggeil l|l992h considered decrease of the 
magnetic field due to Ohmic decay. 

From the statistical ana l ysis o f 24 binary radio pulsars, 
Ivan den Heuvel fc Bitzarakil l)l995l l discovered a clear corre- 
lation between the magnetic field and the mass accreted by 
the NS. In our study, we assume that magnetic field depends 
exponentially on the amount of accre t ed ma tter and use the 
formula suggested bv lOsiowski et all |201ll) : 

-BNS = (-B NS - Bmin) exp(-^^-) + Bmin, (16) 

Mb 

where AM is accreted mass, Mb is magnetic decay mass 
scale and -B m i n is the minimal magnetic field of NS. We 
set M B = 0.025M Q and B min = 10 8 Gauss. B NS is 
initial magnetic field of NS. We assume that the initial 
magnetic fields of nascent NS ar e distributed l og-no rmally 
i|Faucher-Giguere fc KaspH 120061 ; iPopov et all I201QT I. The 
median value and the dispersion of log_B NS are 12 and 1, 
respectively. 

The decay of the magnetic field of a non-accreting NS 
should be s imilar to that of an isolated NS. We assume, after 
iKiel et all l|200ct ). that magnetic field decays exponentially 
due to Ohmic decay 

t tacc \ 



B 



BNsexp(- 



7*B 



(17) 



where t is the age of the NS, tb is magnetic field decay 
time scale and t acc is the time whic h NS has spent a ccret- 
ing matter via Roche lobe overflow. IKiel et all (|2008l ) have 
shown that tb = 2 Gyr is optimal in their preferred models 
of Galactic pulsar population. We use this value of tb in our 
model. 



3 AN EXAMPLE OF NS EVOLUTION 

Symbiotic X-ray bin aries are wide and faint X-ray systems 
ijMasetti et al.ll2006h . In this work, binaries composed of a 
NS and a giant or a giant-like star are considered as SyXBs 
if they satisfy the following conditions: 



(i) The systems are detached; 

(ii) NSs are in the accret or state. 

iRevnivtsev et all (|201ll ) suggested that some LMXBs have 
giant or giant-like donors which fill their Roche lobes, such 
as GX 13+1 and Cyg X-2. Usually, these giant donors are 
still very close to the main sequence and have small He cores. 
In these LMXBs, the X-ray luminosities are very high. Com- 
pared to SyXBs, their orbital periods are very short. In or- 
der to understand SyXBs, we also discuss these LMXBs. 
Further, we use notation "LMXB" for the systems in which 
giant or giant-like donors fill their Roche lobes. 

Below, we present a detailed example of the evolution 
of an accreting magnetized NS in a binary system, as it 
emerges using the algorithm for spin and magnetic field 
evolution presented in the previous sections. We start with 
a binary which consists of a NS and a main-sequence star. 
At the onset of evolution, the mass of NS is 1.4 Mq, its 
surface magnetic field strength is 5 x 10 12 G, and spin period 
is 0.01s. Main-sequence star mass is 1.3 Mq. Initial orbital 
period of the system is 400 day. Figure Q] shows evolution 
of the NS spin period for different input parameters. In 
Table [3] we present some important physical values for all 
points marked by letters in the top left (w w = l/2v asc and 
n = 2) and bottom right (t; w = 2v csc and n — 3/2) panels 
in Fig. [T] For ?; w = l/2iw and n — 2, the evolution of P s 
proceeds as follows. 

Before point A, the secondary is in the first giant branch 
(FGB) stage, but does not fill its Roche lobe. Bondi accre- 
tion rate of the NS is below 1O _15 M0 yr _1 . Neutron star 
evolves like an isolated object. At point A, mass-capture 
rate by NS exceeds 10 -15 Mq yr _1 and NS starts interacting 
with its environment. At point A, relations Rg > Ra > Rc 
(See Table [3} and j a < j'k(-Ra) hold and NS enters the pro- 
peller state of quasi-spherical accretion. While NS evolves 
through the propeller state from point A to point B, its spin 
period increases and _R C increases. At point B, Ra becomes 
lower than _R C , but Mns < 4 x 10 16 g/s. NS enters the quasi- 
spherical subsonic settling accretion state, after which the 
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Figure 1. — Spin evolution of an accreting magnetized NS in a binary system for different input parameters. Light dot-dashed lines 
represent the states in which NS evolves like an isolated NS or it is in the propeller state, that is, the binary does not emit X-rays. Thick 
solid lines represent the states in which the binary is a SyXB. Thick and thin dotted lines correspond to the states in which the binary 
is a LMXB and mass accretion rate onto the NS is lower (between points 'I' and 'J') or higher (between 'J' and 'M') than the Eddington 
critical accretion rate, respectively. In the left two panels, the lines from 'B' to 'E' and from 'H' to 'I' represent the states in which the 
binary is a SyXB. In the right two panels, the binary is a SyXB between points 'B' and 'E'. The inserts in the panels show detailed 
evolution between points 'G' and 'M'. See the text and Table[3]for more details. 



evolution of P a is determined by the torque given by Eqs. 
([TO)) and ©. The stage of SyXB starts. At the beginning, 
i\sd > Ksu, and the spin period is increasing. With the in- 
crease of the spin period, K B a decreases. Soon, at point C, 
A' s d ~ -Ksu, the evolution of the spin period is now deter- 
mined by Eq. (|15p . With the ascend along FGB, stellar- wind 
velocity of the giant w w decreases. This results in the increase 
of j a - At point D, j a > jk(-Ra) and a wind- fed accretion disk 
forms around NS. After point D, the torque exerted on the 
NS is given by Eq {6}. Due to large P s and low Mns, im- 
mediately after point D the spin-up torque is smaller than 
the spin-down torque, and P a increases from 6.3 s to 8.7 s. 
But very soon, with the enhancement of Mns and P a , the 
spin-up torque becomes larger than the spin-down torque, 
and P s starts to decrease. 

At point E the secondary overflows its Roche lobe, 
the system becomes LMXB (in our notation). The mass- 
accretion rate immediately increases and becomes higher 
than the Eddington mass-accretion rate. Neutron star be- 
comes a super-accretor. 

At point F, the secondary leaves FGB and evolves to the 
core helium burning stage. Mass-loss rate rapidly decreases, 
but the velocity of stellar wind is low, so that j a > j'k(^a), 
and neutron star becomes a propeller with a wind-fed ac- 



cretion disk. In this state the mater from the disk does not 
reach the NS due to centrifugal barrier. 

At point G, the secondary evolves to the early asymp- 
totic giant branch (E-AGB). With the ascend along AGB, 
mass-loss rate of the secondary increases. Due to the high 
mass-accretion rate, spin period P s rapidly increases. At 
point H, Ra ~ Rc and j a > j'k(-Ra). The system becomes 
SyXB with a wind-fed accretion disk. At point I the sec- 
ondary overflows its Roche lobe, the system becomes LMXB. 
Very soon, the NS becomes a super-accretor at the point 
J. At point K the secondary evolves into thermally puls- 
ing asymptotic giant branch star. At point M the secondary 
leaves AGB. 

To summarize, the system is a SyXB between points 'B' 
and 'E' and 'H' and T. It is a LMXB from point 'E' to point 
'F' and from point T to point 'M'. 

In our model, in the stages of evolution when the mass- 
accretion rate onto NS is higher than 4 x 10 16 g/s, w w usually 
is low, such that j a > Jk(-Ra) and a disk forms. There- 
fore, the quasi-spherical supersonic accretion hardly occurs 
in SyXBs. 

The case of high wind velocity (« w = 2?; csc ) and respec- 
tive spin period evolution are shown in the right panels of 
Figure [TJ High velocity of the stellar wind from the compan- 
ion to NS results in a low mass-accretion rate by the latter, 
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Table 3. Some important parameters in the turning points of spin evolution marked by the letters in the top left and bottom right 
panels of Fig. [T] The 2nd column gives the evolutionary age of the secondary. Columns 3 and 4 give the masses of the NS and the 
secondary, respectively. The letters in the parentheses in column 4 indicate the evolutionary stage of the secondary: FGB, BHeC, E-AGB 
and AGB stay for the first giant branch (FGB), core helium burning (BHeC), early asymptotic giant branch (E-AGB), thermally-pulsing 
asymptotic giant branch (TP- AGB). The radius of the secondary is given in column 5. Columns 6 and 7 list the terminal stellar wind 
velocity and the mass-accretion rate, respectively. When the secondary overfills Roche lobe, the wind does not play any role. Columns 8, 
9, and 10 give the gravitational capture radius Rq, the magnetospheric Alfven radius Rj± and corotation radius R c , respectively. Column 
11 gives the total torque. Columns 12 and 13 show magnetic field and spin period. 
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J 4945.53 1.44 
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Point Age Mj$s 

(10V) (M ) 

A 4762.05 1.40 

B 4766.22 1.40 

C 4766.23 1.40 

D 4805.38 1.40 

E 4809.80 1.40 

F 4810.18 1.44 

G 4941.68 1.44 

H 4945.39 1.44 

I 4945.47 1.44 
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M 4946.17 1.49 
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which gives a large _Ra- As Figure[T]and Table[3]show for the 
simulation with ii w = 2v csc , at point B of the track when the 
binary becomes SyXB, the age of the optical star is larger 
than the one in the simulation with w w = l/2u csc . In the 
v w = 2u esc case, the binary never becomes a wind-fed SyXB 
when the donor is in the AGB stage, but it may become a 
LMXB for a short time. 

In general, Ra < Rg in the accretor state. The spin-up 
torques in the case of quasi-Keplerian angular momentum 
distribution in the hot shell at the subsonic accretion state 
(ra — 3/2, bottom panels in Fig.[T]), are smaller than those in 
the case of iso-angular momentum distribution (n = 2, top 
panels in Fig. [IJ. Neutron stars in the bottom panels have 
longer spin periods at point C. Therefore, the low value of n 
is favourable to explain SyXBs in which NS have long spin 
periods. 



4 POPULATION OF SYMBIOTIC X-RAY 
BINARIES 

As we mentioned above, theoretical models of the popula- 
tion of LMXBs depend on badly known input parameters, 



such as kick velocity and common envelope treatmen t (e. g., 
Ilben et all Il995l ; iPfahl et alj 120031 ; IZhu et all I20T2T ). How- 
ever, in this pioneering study of SyXBs we focus on the ef- 
fects which are important for the observational appearance 
of the latter: the donor wind velocity (« w ) and the angular 
momentum distribution in the shell of matter settling onto 
NS (index n). 



4.1 The Number and Birthrate of the Galactic 
SyXBs 

About 70% (for «„ = l/2v csc ) — 98% (for « w = 2u esc ) 
of all SyXBs NS are formed via core-collapse, while other 
~ 30% 'K = l/2w csc )— 2% (v w = 2v csc ) via ECS. The 
reason for negligibly low contribution of post-AIC systems 
is large initial separation of components which is necessary 
for formation of an ONeMg WD. The initial orbital periods 
of the progenitors of SyXBs in which an AIC may happen 
are the longest. After formation of the WD, the orbit has 
to remain wide, so that the secondary can evolve to the red 
giant stage (this means that the efficiency of the common 
envelope ejection must be high) . If the orbit is wide, the ac- 
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cretion rate is low and unfavourable for AIC. In the systems 
which produce NS via ECS the primaries usually have ZAMS 
mass between ~ 8 and 11 Mq and short initial orbital pe- 
riods. They encounter a common envelope phase when the 
primaries evolve through HG or FGB. The primaries be- 
come naked helium stars and finally form NS. Core-collapse 
SNe occur in massive stars IIMq) and for progenitors of 
SyXBs with initially massive primaries the range of initial 
orbital periods is limited by the possibility of the system to 
remain bound after SN explosion. 

The input parameter n, describing the angular mo- 
mentum distribution in the hot shell above the NS mag- 
netosphere at the quasi-spherical accretion settling stage, is 
found to have a negligible effect on SyXBs' birthrate and 
number. However, the terminal wind velocity v v greatly in- 
fluences the birthrate and number of SyXBs and introduces 
an uncertainty of up to a factor of about 10. The Galactic 
birthrate and number of SyXB in the simulation with low 
terminal wind velocity, w w = l/2v CBC are ~ 4.1 x 10 -5 yr -1 
and ~ 1000; in the simulation with high terminal wind veloc- 
ity v w = 2v CBC they are lower, ~ 6.6 x 10~ 6 yr~* and ~ 100. 
Up to now, only 10 SyXBs or candidates were discovered. A 
possible reason is that SyXBs are difficult to find because of 
their low X-ray luminosities and transient character of wind 
accretion. 

In contrast to SyXBs, the input parameters n and v w 
are found to have a negligible effect on LMXB (semide- 
tached systems in our notation), population. In the simula- 
tion with v w = 1 /2v CBC , the Galactic birthrate of LMXBs is 
~ 7.0 x 10~ 5 yr -1 , and the number is ~ 10000. In about 50% 
of all simulated LMXBs, NS was formed via core-collapse, 
while in another ~ 50% - via ECS. The contribution of 
post-AIC systems is negligible. Taken at face value, these 
numbers suggest that LMXBs with giant dono rs comprise 
a sign ificant fraction of Galactic LMXBs (cf. IPfahl et al.l 
(2003)). This also means that we encounter the known prob- 
lem of "overproduction" of LMXBs: while population syn- 
thesis codes produce ~ 10 3 — 10 4 strong X -ray systems, th e 
number of observed ones is le ss than 200 (|Liu et alj [2006). 
Ad hoc solution, suggested bv lPfahl et alj (|2003h ispossible 
cyclic behaviour of most LMXB due to irradiation effects, 
thou gh, self-consistent soluti on is not available, as yet. How- 
ever, iBuning fc Ritterl l|2004h found that irradiation does not 
affect giants. The problem remains open. 



4.2 Properties of SyXBs 

Figure presents, in gray-scale, distributions of the orbital 
periods P or b of SyXBs and LMXBs vs. masses of their sec- 
ondary components Mi - Two upper panels show the popula- 
tion of SyXBs, while two lower panels present LMXBs. The 
orbital periods of model SyXBs exceed ~25 days and extend 
to ~30000 days. The range of orbital periods of the known 
SyXBs is from 50.2 day (ICR J16393-4643) to 1161 day 
(GX 1+4). In the distribution for the LMXB population 
there are two peaks of M2 and P rt>. The upper left peak 
is due to the systems with long periods in which the sec- 
ondaries can evolve to the FGB stage and even AGB stage. 
The peak located lower and to the right is formed mainly 
by LMXBs in which the secondaries are in the Hertzsprung 
gap and overflow Roche lobe because of short orbital pe- 
riods. Comparison of the left and the right panels shows 
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Figure 2. — Distributions of the orbital periods P orD of SyXB 
and LMXB vs. masses of their secondary stars. Two upper panels 
are for the SyXBs population, two lower panels for LMXBs popu- 
lation. Gradations of gray-scale correspond to the number density 
of systems >l/2, 1/2 - 1/4, 1/4 - 1/8, 1/8 - of the maximum 
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and blank regions do not contain any stars. 



that high wind velocity is unfavourable for formation of the 
SyXBs and LMXBs with long orbital periods. As Figure [2] 
shows, SyXBs are predominantly the systems with large M2 
and Porb- 

Figure[3]shows the distribution of accretion rates (X-ray 
luminosities) of NS in SyXBs and LMXB. X-ray luminosities 
are approximated as 



L x = n Mc 2 ~ 5.7 x 10 35 erg s" 1 ^) x ( — 



l% B yr-l)> ( 18 ) 



where r\ ~ 0.1 is the efficiency of converting accreted mass 
into X-ray photons. The peak in the distribution of mass- 
accretion rates occurs at ~ 10 _13 Mq yr -1 , and their X- 
ray luminosities are between 10 32 and 10 36 erg s _1 . There- 
fore, SyXBs are faint X-ray sources. Our results are con- 
sistent with observations. LMXBs have high mass-accretion 
rates(~ 10~ 9 M© yr" 1 ). In our work, about 1% of LMXBs 
have super-Eddington accretion. 

Figure [4] shows the distributions of spin periods. The 
range of spin periods in SyXBs is between 0.1 s and ~ 10 5 
s. In the simulation with v„ = 2v CS c and n = 3/2, there are 
two peaks. The left peak is formed by SyXBs with wind- 
fed accretion disks, while the right peak is due to quasi- 
spherically accreting SyXBs because small n is favourable 
for producing long spin periods. In this simulation, about 
20% of SyXBs have wind-fed accretion disks. However, in 
the simulation with w w = l/2v CBC and n = 2, only about 7% 
of SyXBs have wind-fed accretion disks. The low wind ve- 
locity results in a high mass-accretion rate which produces 
shorter spin periods. As the left panel of Fig. Ushows, there 
is only one peak which results mainly from quasi-spherically 
accreting SyXBs. Most of LMXBs have very short spin peri- 
ods (~ 10~ 3 s) because NSs have accreted a large amount of 
matter from their companions. Due to difference in amounts 
of accreted matter, SyXBs have long spin periods, while 
LMXBs have very short spin periods. The latter may be 
good progenitors of binary millisecond pulsars. 

Figure[5]shows the distribution of magnetic fields of NSs 
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Figure 3. — Distribution of accretion rates (X-ray luminosities) 
by NSs in SyXBs and LMXBs. 
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Figure 5. — Similar to Figure [4] but for magnetic fields of NSs 
in SyXBs and LMXBs. 




Figure 4. — Number distribution of the spin periods of NSs in 
SyXBs and LMXBs. The left panel is for population of SyXBs, 
the right panel is for LMXBs. The numbers are normalized to 1. 



in SyXBs and LMXBs. The peak of distribution in SyXBs is 
at 10 11 G. The magnetic fields of NSs in SyXBs are not vio- 
lently suppressed. Most of LMXBs have very weak magnetic 
fields (~ 10 8 G) because NSs have accreted enough matter 



to suppress their fields. 



4.3 Individual SyXBs 



As Table Q] shows, spin periods, orbital periods and X-ray 
luminosities are the most important parameters of SyXBs. 



X-ray luminosities (mass-accretion rates) and spin periods 
of five SyXBs observed and our results for the sample of 
SyXBs and LMXBs are plotted in Figure [6] The models 
with a low wind velocity (« w = l/2« OS c) are preferred for 
explaining the known SyXBs. 

Our results in the simulation with v v — l/2t> eac cover 
the positions of Set X-l, IRXS J180431. 1-273932 and 2XMM 
J174016. 0-290337 very well. Our model cannot identify 
whether there are accretion disks in these systems. Forma- 
tion of the disks greatly depends on the wind velocity. 

Neutron stars in 4U 1954+31 and IGR J16358-4724 
have very long P s . In order to ex plain them by the model 
of a NS with an accretion disk, iPatel et al.l l|2007h sug- 
gested that IGR J16358-4724 contains a magnetar. Indeed, 
as Eqs. 0, (JSJl, and (|15[1 show, the equilibrium NS spin pe- 
riod P s increases with NS magnetic field. However, the po- 
sition of 4U 1954+31 and IGR J16358-4724 in Fig. ©fit well 
the region occupied by quasi-spherically accreting SyXBs 
with low wind velocity. Then an ultra-strong magnetic field 
is not a necessary condition to explain 4U 1954+31 and 
IGR J16358-4724. The same conclusion was obtained by 
iMarcu et all l|201lD who interpreted the 18300 s pulse period 
of 4U 1954+31 within the fra mework of the q uasi-s pherical 
accretion model developed bv lShakura et al.l (|2012l ). 

GX 1+4 has a very interesting long-term spi n 
behaviour dGonz alez-Galan et al. 2012; Shak ura et al.ll2012l ). 
Bas ed on the properties of an assu med accretion disk around 
NS, Ichakrabar tv fc Rochel (|l997l ) suggested that magnetic 
field of NS in GX 1+4 should be ultra-strong (5ns ~ 
10 14 G ). Recently, however, based o n l|Shakura et al.ll2012l ) 
model, Gonza lez-Galan et al.l (|2012| ) argued that GX 1+4 is 
not in the disk-accretion state. Instead, the system which 
currently shows a steady spin-down and can have NS mag- 
netic field close to 10 13 G, definitely is a quasi-spherical 
wind-accreting s ource. As Fig. [6] shows, our re sult supports 
the suggestion of Gon zalez-Galan et al.l l|2012i ). 
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Figure 6. — Gray-scale maps of the distributions of the spin 
period P B vs. mass-accretion rate Mns (or X-ray luminosity L x ) 
for SyXBs. The data on observed SyXBs from Table [l] is plotted 
by cycles. 

As seen in Fig. [6] the main population of Galactic SyXB 
in our population synthesis simulations must have on aver- 
age lower X-ray luminosities than actually observed. This 
difference is most likely due to still very approximate treat- 
ment of the complicated process of stellar wind accretion 
onto NS in our model, which ignores temporal variations of 
the wind properties from cool giants, wind clumping, possi- 
ble effects of the orbital eccentricity, etc. The settling accre- 
tion onto magnet ized NSs at small X -ray luminosities can 
also be unstable l|Postnov et al.ll201l[ ). However, we stress 
good correspondence between observed NS spin periods in 
SyXBs and the ones obtained in quasi-spherical accretion 
model. 



5 CONCLUSIONS 

Using simulation of the interaction of magnetized NSs with 
their environment in binary systems, we investigated Galac- 
tic population of SyXBs. In our simulations, the number 
of Galactic SyXBs is found to be ~ 100—1000, and the 
estimate of their birthrate is between ~ 4.1 x 10 -5 yr _1 
and ~ 6.6 x 10 -6 yr _1 . The simulated SyXBs population 
substantially depends on the properties of the stellar wind 
velocity from cool giants, which is one of the model input 
parameters. In our model, SyXBs have found to have wide 
orbital periods ~(10 — 10000 days), they are faint X-ray 
sources ~ (10 32 — 10 36 erg s _1 ), and have long spin periods 
~(0.1 - 10 5 ) s. 

Our model can explain certain observational properties 
of some known SyXBs or candidate systems. 4U 1954+31 
and IGR J16358-4724, in which NS have very long spin pe- 
riods, are quasi-spherically accreting SyXBs. In this case 
an ultra-strong magnetic field is not necessary condition 
to explain the long spin periods. GX 1+4 also probably is 
a quasi-spherical wind- accreting source. However, we can- 
not identify whether there are wind-fed accretion disks in 
4U 1700+24, Set X-l, IRXS J180431. 1-273932 and 2XMM 
J174016.0-290337. 
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